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ABSTRACT 

We investigate the generation and distribution of high energy electrons in the cosmic structure environ- 
ment and their observational consequences by carrying out the first cosmological simulation that includes 
directly cosmic ray (CR) particles. Starting from cosmological initial conditions, in addition to the gas 
and dark matter related quantities, we follow the evolution of CR electrons (primary and secondary) 
and CR ions along with a passive magnetic field. CR ions and primary electrons are injected in accord 
to the thermal leakage model and accelerated in the test particle limit of diffusive shock acceleration at 
shocks associated with large scale structure formation. Secondary electrons are continuously generated 
through p-p inelastic collisions of the CR ions with the thermal nuclei of the inter-galactic medium. The 
evolution of the CR electrons accounts for spatial transport, adiabatic expansion/compression, and losses 
due to Coulomb collisions, bremsstrahlung, synchrotron and inverse- Compton emission. The magnetic 
field is seeded at shocks according to the Biermann battery model and, thereafter, amplified by shear 
flow and gas compression. 

We compute the emission due to the inverse- Compton scattering of the simulated primary and sec- 
ondary electrons off cosmic microwave background photons and compare it with the published values of 
the detected radiation excesses in the hard X-ray and extreme ultra-violet wave-bands. We find that, 
from the perspective of cosmic shock energy and acceleration efficiency, the few instances of detection of 
hard X-ray radiation excess could be explained in the framework of IC emission from primary electrons in 
clusters characterized by high accretion/merger activity. On the other hand, with the only exception of 
measured flux from the Coma cluster by Bowyer et al. (1999), both primary and secondary CR electrons 
associated with the cosmic structure formation account at most for a small fraction of the radiation 
excess detected in the extreme ultra-violet waveband. 

Next, we calculate the synchrotron emission after normalizing the magnetic field strength so that 
for a Coma-like cluster the volume averaged [B 2 ) 1 / 2 ~ 3/iG. Our results indicate that the synchrotron 
emission from the secondary CR electrons reproduces several general properties observed in radio halos. 
These include the recently found P\aghz vs. T x relationship, the morphology and polarization of the 
emitting region and, to some extent, even the spectral index. In addition, radio synchrotron emission 
from primary electrons turns out large enough to power extended regions of radio emission resembling 
radio relics observed at the outskirts of clusters. Once again we find striking resemblance between the 
general properties of morphology, polarization and spectral index of our synthetic maps and those of 
reported in the literature. 

Subject headings: acceleration of particles — cosmology: large-scale structure of universe — galaxies: 
clusters: general — methods: numerical — radiation mechanism: non-thermal — 
shock waves — X-rays: galaxies: clusters 

1. INTRODUCTION 

The existence of extended regions populated by cosmic ray electrons (CRes) in at least some clusters of galaxies has 
been apparent since the discovery more than thirty years ago of diffuse, non-thermal radio emissions from the Perseus 
and Coma clusters (Leslie & Elsmore 1961; Willson 1970). Their importance as indicators of physical processes in the 
clusters and cluster galaxies has grown in recent years as the number of radio-detected clusters has increased, as reports 
have appeared of possible diffuse non-thermal emissions in the hard X-ray (HXR) and extreme ultra-violet (EUV) bands 
(e.g., Lieu et al. 1996b; Fusco-Femiano et al. 1999) and as the evidence has mounted for a rich variety of highly energetic 
phenomena in and around clusters that seem capable of energizing the electrons (e.g., Sarazin 1999; Enfilin 1999). Early 
on it was recognized that the radiative lifetimes of CRes are almost certainly too short and their diffusion too slow for 
them to fill an entire cluster if they are produced by a single point source (Jaffe 1977). Still, today, the origin of cluster 
CRes is not clear, although many proposals have been made. Our purpose here is to examine through direct numerical 
simulations the possibility that CRes can be explained as a by-product of the very large shocks that accompany the formation 
and evolution of the large scale structure of the Universe. 
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1.1. Radio Observations 

The number of successful detections of diffuse radio emission from clusters of galaxies has substantially increased in 
the recent years and some cluster statistical properties can now be investigated. The radio sources are typically featured 
by low surface brightness and an emissivity that can be described by a steep power law spectrum (Hanish 1982; Enfilin 
et al. 1998). From its spectral properties and, sometimes, polarization signatures the radio emission is interpreted as 
synchrotron radiation, implying the presence of CRcs and magnetic fields (Kim et al. 1989, 1990; Giovannini et al. 1993; 
Deiss et al. 1997; Giovannini & Feretti 2000). Early studies on the origin of the relativistic particles conceived of them as 
originating in discrete radio galaxies (Jaffe 1977; Rcphacli 1977, 1979; Valtaoja 1984) or normal galaxies. More recently, 
the possibility of their acceleration out of the intra-clustcr medium (ICM) thermal electron pool has been considered 
(Liang 1999). Models commonly assume a continual energization of the relativistic electrons by in situ first or second 
order Fermi processes (Jaffe 1977; Roettiger et al. 1993; Tribble 1993; Volk et al. 1996; Deiss et al. 1997; Eilek & Wetherall 
1999; Wetherall & Eilek 1999). Alternatively, such electrons could be produced as secondary products of inelastic p-p 
collisions of CR ions and the thermal intra-cluster nuclei (Dennison 1980; Vestrand 1982; Blasi & Colafrancesco 1999). 
This idea was put forward primarily to explain the wide-spread distribution observed for these electrons, which is in 
conflict with the fact that their cooling time scale is too short for them to propagate throughout a cluster from a single 
or even a few discrete sources (Dennison 1980). 

Depending on its observed properties the diffuse radio emission is further classified as either a radio halo when its 
morphology is regular and typically centered on and resembling the X-ray emissivity, or as a radio relic when it is 
irregular and located at the periphery of the cluster. The former typically has a low degree of polarization (< 10%). The 
degree of polarization in the latter case, when measured, turns out typically to be high (~ 30%; e.g., Feretti 1999). 

Radio halos are usually found in rich clusters with high ICM temperature, TJ>7keV, and high X-ray luminosity, L^O.l — 
2.4keV)>5 x 10 44 erg s _1 (Feretti & Giovannini 1996). Since it usually extends over a linear size of about 1 /i _1 Mpc, the 
radio emission appears to be a characteristic of the whole cluster, rather than being associated with any of the individual 
cluster galaxies (Willson 1970). Signatures of a merging process in these clusters are often emphasized and the absence 
of cooling flows cited as demonstrating the connection between radio halos and mergers (a merger would likely disrupt 
a cooling flow; e.g., Feretti 1999). This idea seems attractive, because a merger can provide enough energy and the 
turbulence necessary for the amplification of the magnetic field and possibly acceleration of relativistic particles needed 
to power a radio halo (Tribble 1993). In addition, since mergers are transient events, they would also explain why this 
type of radio source has appeared to be only occasional. Recently, however, Liang et al. (2000) found a tight and steep 
correlation between the radio power emitted at 1.4 GHz and the cluster temperature, suggesting that the apparent rarity 
of detections should be attributed to observational insensitivity to any but the rarest, most massive clusters. Those 
authors also claim that the radio properties of clusters with cooling flows have not been explored sufficiently as yet. 

On the other hand the clusters hosting radio relic sources are somewhat less massive and cooler than those related to 
radio halos. Radio relics show no apparent correlation with merger events (although the two can be coincident; cf. Feretti 
1999) and they are definitely observed in clusters containing cooling flows (Bagchi et al. 1998). Radio relics are found 
both near the cores of clusters and at their outskirts. The spectra are typically steep, but explicit cutoffs are relatively 
rare even though the cooling time of the relativistic electrons is much shorter than the age of the cluster. Thus an efficient 
production mechanism for an extended particle population is required in order to explain the origin of those relativistic 
electrons, as well. In Enfilin's model (Enfilin 1999) a seed population of relativistic electrons is provided by radio ghosts, 
i.e., relic relativistic plasma previously injected by radio galaxies. Such a plasma would be re-energized by the encounter 
with a shock due either to a merger event or to supersonic inflow along cosmic filaments (Enfilin et al. 1998). Given the 
high sound speed inside the relic plasma, such flows would actually become sub-sonic and propagate as a sound wave 
through the relic, mostly producing only an adiabatic compression of it (Enfilin & Gopal-Krisna 2001). 

1.2. Hard X-ray and Extreme Ultra-violet Observations 

An energy excess with respect to the thermal emission has been detected in the HXR band of the spectrum at least 
in four instances: for the Coma cluster (Fusco-Femiano et al. 1999), for Abell 2199 (Kaastra et al. 1999) and for Abell 
2256 (Fusco-Femiano et al. 2000). An upper limit for a HXR non-thermal excess was recently measured also for Abell 
3667 (Fusco-Femiano et al. 2001). Although possibly originating from high temperature shocked intra-cluster gas, HXRs 
are also expected as a consequence of inverse- Compton (IC) scattering off cosmic microwave background photons by 
electrons responsible for synchrotron emission. An approximate "RMS" magnitude of the magnetic field can be inferred 
from combining HXR and synchrotron surface brightness estimates, under the assumption of a common population of 
relativistic particles generating the two observed radiations. However, the result for the Coma cluster indicates a value of 
B <~ 0.16/iG (Fusco-Femiano et al. 1999), which is more than an order of magnitude lower than magnetic field estimates 
through Faraday rotation measures, B <~ 6/iG (Feretti et al. 1995) (see also Clarke et al. 2001). An equipartition estimate 
of magnetic field and relativistic electrons results in B ~ 0.4^iG (Giovannini et al. 1993). There is no good theoretical 
argument for the existence of that kind of equipartition, however. Given these uncertainties, the origin of the HXRs is 
still open to debate. For this reason, non-thermal brcmsstrahlung emission has also been proposed to explain the origin of 
the HXRs (Enfilin et al. 1999; Sarazin & Kempner 2000; Blasi 2000; Dogiel 2000). However, Petrosian (2001) argues that, 
because of the very low efficiency of non-thermal bremsstrahlung with respect to Coulomb losses, in the above model a 
large amount of energy would need to be dumped in the ICM, leading to an unacceptably large consequent equilibrium 
temperature. 
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Additional evidence for the existence of non-thermal populations of rclativistic electrons may be suggested by obser- 
vations of clusters in the EUV. There have been reported detections of EUV radiation from clusters in excess of what 
is expected from the hot, thermal X-ray emitting ICM (Lieu et al. 1996b, a; Fabian 1996; Mittaz et al. 1998; Lieu et al. 
1999; Kaastra 2000; Bonamente et al. 2001). This result appears, however, to be sensitive to the procedure adopted for 
the subtraction of the background (see, e.g., Lieu et al. 1996a; Mittaz et al. 1998; Bowyer et al. 1999, for the details). 
In any case, one mechanism proposed for the origin of an EUV excess is IC scattering off cosmic microwave background 
photons of low energy relativistic electrons with a Lorentz factor of 7 <~ 300 (Sarazin & Lieu 1998). As an alternative, the 
possibility that the EUV originates from warm gas produced at the interface between a population of cold clouds and the 
hot ICM has also been explored (Fabian 1996; Bonamente et al. 2001). However, observations of the FUV spectra with 
FUSE has ruled out the existence of a substantial amount of warm gas, at least, in the Coma and Virgo clusters (Dixon 
et al. 2001). 

In summary, there are indications of excess, "non-thermal" HXR emissions and, possibly EUV emissions, from some 
clusters that may signal the presence of substantial CRe populations there. This information compliments the radio data, 
even in its current, uncertain state, and adds to the motivations for understanding the likely sources and distributions of 
CRe populations in clusters. 

1.3. Cosmic Shock Waves 

Collisionless plasma shocks are now generally recognized as efficient CR accelerators through the so-called "diffusive 
shock acceleration" (DSA) process (e.g., Drury 1983). Cosmic structure formation simulations have demonstrated that 
strong accretion shocks form several Mpc's from cluster cores and also penetrate deeply into cluster cores during the course 
of large-scale structure formation (Kang et al. 1994; Miniati et al. 2000). They have been suggested as possible sources of 
high energy CRs (Kang et al. 1996, 1997) as well as seeds for the ICM magnetic field (Kulsrud et al. 1997; Ryu et al. 1998). 
Miniati et al. (2000, 2001) pointed out that those shocks associated with large scale structure formation, in principle have 
sufficient power to produce dynamically significant CR proton populations inside clusters, and by implication substantial 
populations of CRes. Those computations, based on cosmological simulations, also showed, however, that cluster shock 
structures are very complex. Thus, it is necessary to explore in some detail the physics of CR acceleration in this context 
before one can properly evaluate its role. In this paper we will investigate by means of a numerical simulation CRes 
injection and acceleration by such cosmic shock waves. In particular, we explore the possibility that the origin of the 
aforementioned radio emission and radiation excesses can be explained by the primary and secondary CRes and focus our 
attention on those quantities that are directly related to observations. These include the flux and the spectral features 
of the radiation emitted at various wavelengths, including the EUV, soft and HXRs produced by IC scattering and radio 
synchrotron emission from the CRes. 

There have been previous attempts to study the observational consequences of CRes with simulations (e.g., Roettigcr 
et al. 1999a; Takizawa & Naito 2000; Dolag & Enfilin 2000). However, for the first time we include the CR population 
explicitly in a fully cosmological simulation by computing particle injection, acceleration and energy losses in accord with 
the properties of the local environment in which the particles are propagating and starting from cosmological initial 
conditions. Here, our focus is CR "electrons"; CR "protons" (as the dominant ionic CR component) are discussed in a 
companion paper (Miniati et al. 2001). We note the suggestions that there are additional sources of CRs in clusters, of 
course, such as active galaxies (Volk et al. 1996; Enfilin et al. 1997; Berezinsky et al. 1997). We do not attempt to include 
them in our current work, since our goal is to understand the role of structure shocks. However, we do call attention in 
our discussions to some expected differences between shock CR sources and isolated sources as appropriate. Our modeling 
efforts provide some initial clues as to how these different sources can be distinguished observationally. 

The paper is organized as follows. In section 2 we describe the numerical methods. In §3 we present the main results 
of our simulation and the conclusions are given in §4. 

2. THE NUMERICAL SIMULATION 

In order to investigate the various issues related to CRes, we have carried out a cosmological numerical simulation that 
follows simultaneously the evolution of dark "matter" particles, gas quantities, magnetic field, plus CR ions, primary 
CRes injected from thermal plasma and secondary CRes. The cosmological aspects of the simulation setup are presented 
in the following subsection (§2.1). For the matter part, we have used an Eulerian "TVD" hydro + N-body cosmology 
code (Ryu et al. 1993). The magnetic field is seeded at shocks in accord to the Biermann battery mechanism and is 
thereafter followed passively (Kulsrud et al. 1997). After being generated at shocks, the magnetic field is amplified by 
shear flow through stretching and by field compression. Numerical resolution in this simulation, however, is insufficient 
to produce field strengths of the same order of magnitude as observed in today's groups/clusters of galaxies from realistic 
seed fields. Therefore, the strength of magnetic field is normalized a posteriori to be consistent with the values inferred 
from observations of Faraday rotation measure through clusters (e.g., Clarke et al. 2001). Finally, the evolution of the 
different CR populations in the simulation is computed by the code COSMOCR (Miniati 2001). In the following sections 
we provide a brief description of the physical processes included in this code, i.e., CR injection at shocks and spatial 
transport (§2.2), production of secondary electrons (§2.3) and energy losses (§2.4). 
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Since the computation involves a quantity never before simulated in this context, i.e., the CRs, we have elected to begin 
the study from the relatively simpler case of a standard cold dark matter (SCDM) model, based a modest resolution 
simulation with 256 3 cells and 128 3 dark matter (DM) particles. We leave the currently more favored CDM + A model as 
the natural follow-up step for future work. Although it is now well known that SCDM is not a viable model (e.g., White 
et al. 1993; Ostrikcr 1993), the key feature required of the simulation is a distribution of collapsed objects whose properties 
resemble those observed in the real universe, thus allowing comparisons of their general characteristics. To accomplish 
this within SCDM, we set the value for the root mean square of today's density fluctuation on a scale of 8 h~ 1 Mpc to be 
(7s = 0.6, incompatible with COBE results and a SCDM universe, yet a value that induces the emergence of a reasonable 
population of groups/clusters in simulations of large scale structure formation (Ostriker & Cen 1996). In addition to cr 8 
we also define the following parameters: the spectral index for the initial power spectrum of perturbations, n = 1, the 
normalized Hubble constant, h = H o /(100 km s _1 Mpc -1 ) = 0.5, the total mass density Qm = 1 and the baryonic mass 
fraction fif, = 0.13. 

The size of our simulated region is chosen as a compromise between a simulated cosmological volume large enough to 
contain groups/clusters of galaxies and a resolution sufficient to capture the main structural features of the simulated 
objects (groups/clusters and their shocks). We have selected a cubic region of size 50 /i _1 Mpc at the current epoch. 
With 256 3 cells this corresponds to a spatial resolution of about 200/i~ 1 kpc. In general a coarse grid limits the finest 
structure that can form and consequently influences properties of the simulated objects. Then, because density peaks 
are smoothed, estimates of any quantity depending on the square of the density, such as the thermal brcmsstrahlung 
emission (oc n^ hermal ), will be reduced. The effect is stronger for lower temperature groups/clusters, which have similar 
structures to the larger clusters, but also have smaller physical scales. Thus, because of resolution effects the emissivity 
from secondary electrons (whose production roughly scales as oc nf hermal ) will be systematically underestimated and will 
have a steeper temperature dependence. Previous resolution studies carried out to test the performances of the code 
employed here indicate that the thermal emission is underestimated by a factor of a few, although the order of magnitude 
should be correct (Cen & Ostriker 1999). On the other hand, the resolution errors on the emissivity from primary electrons 
should be less severe compared to the case of secondary electrons. 

2.2. Injection at Shocks and Spatial Transport 

Ions and electrons injected at shocks from thermal plasma should provide the main source of CRs and are usually referred 
to as primary CRs. However, the hadronic interaction of high energy ions with the thermal nuclei of the background 
gas generates, among other secondary decay products, relativistic electrons and positrons. The COSMOCR code follows 
explicitly the injection, acceleration, transport and energy losses of three CR species: ions, primary electrons and secondary 
electrons. The complete CR computational methods are described in Miniati (2001), and the detailed results of the CR 
ion simulation are presented in Miniati et al. (2001). 

The injection of CR ions at shocks is computed in the "test particle limit" and according to a "thermal leakage" 
model (Kang & Jones 1995). The fraction of supra-thermal protons injected into CRs from plasma flowing through 
a shock is c ontrolled by an injection parameter, a, which defines the momentum threshold for injected protons as 
Pinj — c i 2^/TOpfcTs/j. In this expression m p is the proton mass, ks the Boltzmann's constant and T s h the post-shock 
plasma temperature. We adopted c\ = 2.6, which leads to injection of <~ 10~ 4 of the thermal protons passing through 
shocks. This injection rate is consistent with observation and theoretical studies of DSA (Lee 1982; Quest 1988; Kang & 
Jones 1995; Malkov 1998). 

Since thermal electrons have gyro-radii much smaller than analogous ionic gyro-radii, they cannot be injected in the 
same fashion as protons. A pre-acceleration process is needed before electrons can be introduced into the DSA process 
(e.g., McClements et al. 1997). However, the physics underlying the injection at shocks of primary electrons is very 
complex and not well understood. Therefore, following the approach of Ellison et al. (2000), we compute the injection 
rate of CRes by assuming a fixed value R e / P for the ratio of primary electrons to protons at relativistic energies. There is 
some observational evidence that this ratio is in the range 1-5 % for the Galactic CRs (Mueller & Tang 1987; Mueller & 
et al. 1995). Thus, the normalization of the number density in our calculations is controlled through the ionic injection 
rate for secondary electrons and both the ionic injection rate and R e / p for primary electrons. For convenience we report 
our results with R e / P — 1. Thus the emissivity for the primary electrons needs to be multiplied by an additional factor 
Re/p ~ 0.01 — 0.05 in order to get practical values. 

Spatial diffusion of the GeV CRes on the scales of interest here occurs on time-scales much longer than advection and 
can therefore be neglected in regions dominated by smooth flows, i.e., shock-less regions (Jones et al. 1999; Miniati 2001). 
However, the diffusive properties of the CRs are still important because: (1) the acceleration process depends explicitly on 
diffusive propagation immediately adjacent to the shock and both (2) the maximum energy reachable by the accelerated 
particles and (3) the maximum energy of the particles confined within the ICM depend on the diffusion coefficient. As for 
point (1), except for extremely high energy CRs the scales on which that diffusion applies are far too small to be modeled 
in a cosmological simulation, and the time scale for the acceleration is effectively instantaneous by the same measures. 
Therefore, as explained in Jones et al. (1999), the injected CR protons and electrons are redistributed in energy according 
to a power law as prescribed by DSA theory (e.g., Drury 1983). Analogously, a pre-existing population of CRes passing 
through a shock is compressed, and, if appropriate, its spectrum flattened according to the Mach number of the shock. 
Additional details of these procedures can be found in Miniati (2001). As for points (2) and (3), the diffusion coefficient 
and the associated mean free path are usually computed in the pitch angle resonant scattering approximation (e.g., Drury 
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1983) as 

1 B 2 /8tt 

where W\- is the magnetic energy density per wavenumber evaluated at k <~ l/r Ll i.e., the inverse of the Larmor radius. 
Thus the diffusion coefficient is inversely proportional to the level of resonant wave turbulence, and depends, therefore, 
on the spectral properties of the turbulent magnetic field. The diffusion coefficient determines the rate of energy gain 
which, when compared with the dominant energy loss rate gives an estimate of the maximum energy of the accelerated 
particles. It is shown that in the case of an average ICM magnetic field of value <~ 0.15^G and with a Kolmogorov power 
spectrum, the upper limit for the energy of the accelerated electrons could be low enough to produce a cut-off in the radio 
synchrotron emission around a few GHz (Blasi 2001). Such a cut-off, however, disappears if one makes the assumption of 
Bohm diffusion (Blasi 2001). In any case, with the normalization adopted in §2.4 the magnetic field in the groups/clusters 
of interest here (with temperature T x > IkeV) is always larger than the above value of 0.15 ^G. Therefore, such a cut-off 
should not affect the spectra of our simulated clusters. Finally, we point out that diffusive escape of the CR-ion component, 
which is responsible for the production of the secondary electrons, only affects particles with energy above ~ 2 x 10 15 eV 
(Volk et al. 1996; Berezinsky et al. 1997). This component was not included in our simulation (Miniati et al. 2001) and 
should not be relevant for the secondary electrons in the energy range investigated here (e.g., Mannheim & Schlickeiser 
1994). 

2.3. Production of Secondary Electrons 
The main channels for the production of secondary electrons are given by (Gaisser 1990): 

p+ P ->TT ± +X (2) 

p + p^K^ + X, (3) 

where X indicates all the other by-products of the reactions. Subsequently charged pions and kaons decay mainly into 
muons (although a fraction of kaons also produce pions), from which the secondary electrons originate (see Miniati 2001, 
for more details on the processes included in COSMOCR). In addition to p + p inelastic collisions, the above cascades are 
also triggered by the interaction of p+He, a+H and a+He which, for example, increase the the overall yield of secondary 
e~ by a factor 1.4 for a metal composition similar to the inter-stellar medium (Dermer 1986). Here we have assumed a 
helium number fraction of 7.3% for the background gas and a ratio (H/He) ~ 15 at fixed energy-per-nucleon for the CRs 
(Meyer et al. 1997). 

In general we write the production spectrum of secondary electrons as (Moskalcnko & Strong 1998): 

/ ds p J p {s p ) (Oi(e P )) / d£iF a {e s ,£i)Fi{£i,£p), (4) 



=tt,K- 



where J P (£ P ) is the proton flux; ((ai(£ p )) is the inclusive cross section 7 of the processes 2 and 3; £™ m is the minimum 
proton energy required to produce a meson of energy e™ m and £™ aa: the maximal energy of the produced meson; £™ nn in 
turn is the minimum required energy of a meson for production of secondaries of energy e s ; finally Fi(e ni K, £ p ) are the 
spectra of n and K produced from the collision of a proton of energy £ p and F s (£ Sl £ 7r ^) the distribution of secondaries 
from the subsequent decay of the above collision products. The cross sections were computed according to the GALPROP 
routines by Moskalcnko & Strong (1998) to which we refer for a full summary of the technique employed here for the 
calculation of the secondary electrons. 

2.4. Radiative Losses 

COSMOCR (Miniati 2001) accounts for various energy loss mechanisms suffered by the CRcs in different energy ranges 
of interest here. To do this, a Fokker-Planck equation that has been integrated over finite momentum bins is solved. This 
takes advantage of the near-power-law form of the CR momentum distribution. In effect, the momentum space is divided 
into N p = 8 logarithmically equidistant intervals, bounded by p , ...pN p = Pmax, which we refer to here as momentum 
bins. The electron distribution function f(p), as a function of the normalized momentum p = p/m e c, in each spatial cell 
and for each momentum bin is approximated by the following piece- wise power law: 

/(Xj.p) = fj(xi)p-«™, 1 < Pj-i <P< Pj, (5) 

where /j(xj) and (?j(xj) are the number normalization and logarithmic slope for a given cell and p bin. With the above 
definition, the number density of particles is given by dN — Airp 2 f(p)dp. 

Within each momentum bin, j, and at each spatial grid point, x i; we follow the CRes total number density and kinetic 
energy density defined as 

/•Pj+i 

n(x i ,p j )=4Tr / f(xi,p)p 2 dp (6) 
Jpj 
rPi+i 

E(xi,pj)=4ir / ,f(x u p)T(p)p 2 dp, (7) 
Jpj 

7 Inclusive means that which describes the process p + p — > i + X, where i is in general a secondary particle. 
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where T(p) = (7 — l)m e c 2 is the relativistic kinetic energy. Further, for each momentum bin qj(xi) is determined self- 
consistently from the values of n(x i ,p J ) and s(xi,pj) defined above (see Miniati 2001, for details). With this formalism, 
the evolution of n(x i ,p J ) in momentum space is described by the equation 

9n{X g f j) = -V • [unfxi.ft)] + [b(p) 47TP 2 m^+Qi^pj) (8) 

where the first term on the right hand side describes advective transport and Q(pci,pj) represents the source term, i(xj,j5), 
integrated in the jth bin. Finally, b(p) = dp/dt describes mechanical and radiative loss terms (see Miniati 2001, for further 
details on this). The most effective among these are Coulomb losses in the low energy end and synchrotron and IC 
emission at the highest energy. Bremsstrahlung losses are less important, although included for completeness. Losses due 
to Coulomb collisions are (Strong & Moskalcnko 1998) 

(ij^-uxio-ji + Ki+^-M^j}-^. 

where n is the number density of the background gas in cm~ 3 . Bremsstrahlung losses are defined as (Ginzburg 1979; 
Strong & Moskalenko 1998) 



dt 



1.4 x 10~ 16 |ln[2 (1 + p 2 ) 1/2 } - ^ j pn s _1 . 



brem 

For synchrotron emission and IC scattering, the combined contribution is given by 



(f) = 1.36 x 10" 20 (l + ^ W\ 

V dt J svnc+IC \ Ucmb ) 



sync+IC 

where Ms is the energy density in magnetic field, and u cm b = 4.2 x 10~ 13 (1 + z) 4 erg cm -3 is the energy density in cosmic 
microwave background. 

Since our simulation treats the magnetic field passively, it is necessary to normalize the magnetic field strength. For 
this purpose, we assume that the volume averaged magnetic field in a Coma- like cluster of <~ 8 keV is of order 3 /iG, 
which corresponds to values inferred from Faraday rotation measures (Kim et al. 1989, 1990). Thus for most of the 
groups/clusters in our simulation, characterized by a temperature T x < 3keV, the energy density in the magnetic field is 
smaller than that in the cosmic microwave background and the radiative loses due to IC process dominate over synchrotron 
losses. The IC cooling time of CRes is of order 

T ic = A.i> X 1U I ■ 

where p — 7 ~ 300 (10 4 ) is required for electrons generating IC emission in the EUV (HXR) band. 



t; c = 2.3x10 8 (^) yr, (9) 



2.5. Calculation of Emissivity 

In this section we described the procedure adopted for the calculation of IC emission in the HXR and EUV bands as 
well as radio synchrotron at various wavelengths of interest. 

2.5.1. Inverse- Compton Emission 

IC emission is computed in units of erg cm~ 3 s _1 Hz -1 sr _1 according to (e.g., Rybicki & Lightman 1979, but note 
the different units) 



3 cut 



N p 



3ic{v) = " Y 1 hv I v cmb (v )d\nv ^fj f g ( ^ \ p q] dp (10) 
4 J j—i Jpj-i V Q J 

where g{x) = 2x \nx—2x 2 +x+l, w cm 6(fo) = u cmb {v )/hv indicates the number density of photons in the radiation field at 
frequency v and the sum extends over all the bins of the simulated electron populations defined in 5. The above equation 
10 is the actual formula that we used to calculate the inverse Compton emission at various wavelengths from the simulated 
electrons. Since most of the emission is contributed from the CRes in the k t h bin such that Pk-i < (^/4^o) 1//2 < pk, for 
the sake of illustration we also write (Rybicki & Lightman 1979) 

j ic {u) * f*-^pt J L(q k - 2) ,-(-- 3 )/ 2 . (11) 

where 

L (a) - 2^+ 3 g2 + 4g+U r ( ^] V n-(«+ 5 )/ 2 f 12) 

{q) ~ ( (? + 3) 2 ( (? + 5)( (? + l)U 2 ■ (12) 

Thus from cq. 11 we infer that the frequency dependence of the IC emissivity is a power-law, j(v) oc v~ a , with spectral 
index 

a~— — . (13) 
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2.5.2. Synchrotron Radiation 

Analogously, for a piece- wise power-law distribution of electrons given in eq. 5, the synchrotron emissivity (in units of 
erg cm~ 3 s _1 Hz -1 sr _1 ) is computed according to the expression 

where = eB / (2nm e c) , Xj = \.hv /vbPJ 2 and we approximate F(x) = 0.87a; 1 / 3 exp (-lla; 7 / 8 /8) (Enfilin ct al. 1999). 
As for the inverse Compton emission we can also derive an approximate expression (Ginzburg & Syrovatskii 1965) 



K A 4^ r (» + H V (it _ n c^y^- 3 " 2 _ (15) 



with k referring to the momentum bin such that f>k-i < y/^.hv/vB < Pk- Again, from 15 we infer a spectral index 
a ~ (qfe — 3)/2, as for the IC emission case. The local magnetic field strength, entering vg, is given by the simulated value 
after a global re-normalization is performed according to the criterion given in §2.4. Further, by using the morphology of 
the simulated field, the degree of polarization of the synchrotron radiation is also computed. 

2.6. Data Analysis 

From the simulation results collapsed objects have been identified by the DM-based "spherical over-density" method 
described in Lacey & Cole (1994). The details of our identification procedure can be found in Miniati et al. (2000). 
Global quantities, Q, (e.g., luminosity, photon flux, spectral index), are extracted by integrating or averaging over the 
group/cluster volume as 

Q = ^2 Qi w i ( 16 ) 

i 

where qi is the cell value of the quantity of interest and Wi is the weight function, i.e., the portion of computational 
cell within V = Values of Wi have been computed up to a diameter of 13 computational cell corresponding, 

given our current resolution, to radius R c i ~ 1.3 h~ 1 Mpc. Some of these global quantities are studied as a function 
of the group/cluster core temperature, T x . The latter is defined as the average temperature within a central region of 
diameter 0.5 /i -1 Mpc, the typical size of the observed X-ray clusters (and over which the temperature is observationally 
determined). A temperature dependence is evaluated in general by fitting a computed global quantity, Q, to the following 
power-law, 

Q = K ( Tx \ . (17) 
* V6.72kcVy v ; 

where K and <j> are the fitting parameters. The best fit parameters are then retrieved by a least x 2 analysis. Global 
quantities such as the luminosity at different wavelengths are expected to correlate with T x since the emitted radiation 
should depend on the available group/cluster gravitational energy. In any case, we use the values of the scaling parameters 
to make predictions about the quantities under investigation and to make comparisons between the predictions from our 
simulation and observations. 

Alternatively, two-dimensional maps have been constructed from the data-set, either simply as slices through the volume 
of physical a physical quantity (e.g., density); or by way of a line of sight integration of a volume emissivity relative to 
a selected simulated group/cluster. The latter procedure produces synthetic surface brightness maps in the optically 
thin approximation through a projection code (Tregillis 2001). For comparison purposes the groups/clusters are set to 
a luminosity distance corresponding to the red-shift of the Coma cluster (z — 0.023), i.e., about 70 /i Mpc. Since our 
resolution amounts to ~ 200 /i _1 kpc per computational cell, at this distance the minimal size of a pixel of the synthetic 
image corresponds to 9.8 arcmin squared. These images allow a more in-depth inspection of the spatial distribution of 
the quantities of interest, but for practical reasons are limited so far to only a few samples. 



3. RESULTS 
3.1. Inverse- Compton 

IC emission from the simulated CRes is computed in the following wave-bands: in HXR, between 20 and 80 keV, where 
measurements have been carried out for Coma cluster (Fusco-Femiano et al. 1999), Abell 2256 (Fusco-Femiano et al. 2000) 
and Abell 3667 (Fusco-Femiano et al. 2001), and between 0.13 to 100 keV, relevant for Abell 2199 (Kaastra et al. 1999); 
and in the EUV between 65-248 eV, the wave-band of the EUV Explorer satellite. 

Surface brightness maps for one of the collapsed objects generated in the simulation are presented in Fig. 1 for the HXR 
emission between 20 and 80 keV in units 2 x 10~ 4 erg cm~ 2 s _1 per pixel and in Fig. 2 for the EUV emission in units 
4 x 10~ 7 erg cm~ 2 s _1 per pixel, respectively. In each figure, left and right panels correspond to secondary and primary 
electrons, respectively. From Figs. 1 and 2 it appears that the radiation from primary electrons is more extended, less 
centrally peaked and morphologically more irregular than that from secondary electrons. The irregularity of the maps 
corresponding to the primary electrons is caused by the fact that latter are preferentially found around existing shocks, 
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where they are freshly injected and have not yet been depleted by radiative losses. Such irregularity is stronger for HXR 
than for EUV photon energies. In fact, primary HXR emitting electrons, with Lorentz factors 7 <~ 10 4 (cf. eq. 9), have 
very short cooling times and do not propagate far from the shock where they have been generated. For the primary EUV 
emitting electrons, the cooling time is about 100 times longer, and in fact these electrons are distributed throughout the 
collapsed object, although they arc still more numerous around shocks. The distribution of secondary electrons, on the 
other hand, and the ensuing morphology of IC emission maps, are much less affected by, but not totally insensitive to, the 
presence of shocks. Their emissivity roughly varies as n 2 gasl indicating that although spread more widely, the distribution 
of the parent CR protons follows more or less that of gas (Miniati et al. 2001). Such general properties of the emission 
from primary and secondary electrons will characterize the radio maps as well. 

The total IC emission in the wave-bands of interest has been computed for each simulated group/cluster and plotted 
as a function of its core temperature, T x . In Fig. 3 we report the fluxes of HXR radiation in the ranges 0.13-100 keV 
(left panels) and 20-80 keV (right panels) produced by secondary (top) and primary (bottom) electrons respectively, for a 
luminosity distance of the simulated groups/clusters of 70 /i _1 Mpc. In Figs. 4 we present the luminosity in the EUV band 
(65-248 eV) for secondary (left) and primary (right) electrons respectively. The values reported for the primary electrons 
assume R e / P = 1 (sec §2.2 for further detail). The data plots have been fitted to a power-law correlation as described in 
§2.6. The best fit parameters are given in Table 1. Note that with primary electrons the slope <$> for 20-80 keV HXR is 
somewhat smaller than those for 0.13-100 keV HXR and EUV. This is at least partly because of larger scatter seen in the 
bottom-right panel of Fig. 3. We believe that this scatter has a physical origin in the different shock acceleration history 
undergone by each collapsed object identified in the simulation box. However, we do expect larger simulations, which 
allow both more numerous and larger clusters as well as higher resolution, to result in more consistent values of <f>. 

The temperature scaling of the emitted radiation from the two CRe populations is quite different. This fact corresponds 
to differences in the scaling properties of secondary and primary electrons with the group/cluster temperature. The scaling 
of the IC flux from secondaries is expected to go as 

F ic (v) oc N es (^ v ) oc N cr (^ v ) n gas , (18) 

where iV e s (7i/) an d N cr (-f u ) are the total number of secondary electrons with energy 7^ responsible for the emission at the 
specified frequency, u, and that of the parent CR ions respectively. By accounting for the gas accretion rate and for the 
injection mechanism adopted here as roughly converting a fixed fraction of the shocked gas into CRs (see §2.2), it can be 
shown that N cr {^ u ) ccT^. In addition, the average baryon density inside the group/cluster volume within the fixed radius 
of 1.3 /i _1 Mpc has an almost linear dependence on T x , in accord with observations (Edge & Steward 1991; Mohr et al. 
1999). Although resolution effects discussed in §2.1 may cause some steepening of the Fi C — T x relation, our numerical 
results are consistent with the above arguments. On the other hand, for the IC emission from primary electrons we have 

F tc {v) cx N ep { lv ) ~ N cr { lv ), (19) 

since the total number of primary electrons and ions should be proportional at the time of injection. Therefore, in this 
case the temperature dependence of the flux from IC scattering should be weaker by roughly one power of T x with respect 
to the secondary CRes model, as indeed we have found. 

We have also studied the spectral properties of the radiation emitted by the simulated electrons. Since the published 
data report observed values of the spectral index for the HXR emission only, we only present our results for this wave-band, 
although the same general discussion described below applies to EUV emitting electrons too. Thus, for each simulated 
group/cluster we have computed the IC HXR spectral index, oti C , averaged over a cluster-centered volume defined by a 
radius R to give cti C {< R). Our results are reported in Fig. 5 for both primary (bottom panels) and secondary electrons 
(top panels), as a function of the clusters/groups temperature, T x , for four different values of R. 

As for the primary electrons the situation is complicated because our spatial resolution is inadequate to follow in detail 
the spectral evolution of this population of CRes immediately downstream of shocks. In fact, after being shock-accelerated, 
CRes will rapidly radiate away their energy as they are advected away from shock acceleration regions. Thus, in a realistic 
picture, at various distances from the shock we should find populations of CRes with increasingly "older age" . However, 
for initial energies of interest the time-scale for cooling is typically shorter compared to the time for advection across 
a numerical cell. Therefore, the proper CRes spatial distribution cannot be reproduced because of resolution limits. 
Nevertheless we know that for the case of IC and synchrotron cooling, the summation of all of these "aging" populations 
in the post-shock region produces a volume- averaged population with a power-law index steeper by one from that produced 
at the shock. Thus, the spectral index of the volume-averaged emission from the post-shock region can be recovered from 
the properties of the shock. In addition, in order to compute correctly the total emission from the post-shock region, we 
adopt a sub-cycling procedure that allows us to estimate the steady-state number density of high energy electrons there. 

With this clarification, in the four bottom panels of Fig. 5 we present the volume-averaged spectral index associated 
with the post-shock IC emission region. According to our discussion above, the spectral index characterizing the IC 
emission from the high energy electrons that are emerging from the acceleration region, i.e., those immediately behind 
the shock is flatter by about 0.5 (cf. eq. 13). For shocks with Mach number > 3, representing most of the shocks in 
our simulation, the spectral index is oti C < 1.3. Indeed, in the last panel (bottom-right) of Fig. 5, corresponding to an 
averaging volume of radius 1.3 /i Mpc around the group/cluster center we have cti C ~ 1. — 1.3. This volume is large 
enough to include most of the shocks where fresh particles are injected. On the other hand, for smaller radii, most of the 
considered volume has been depleted of energetic particles through radiative losses. Therefore, we find higher values and 
larger fluctuations for a, c (< R). 
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Turning now our attention to the top panels of Fig. 5, we notice that for the secondary electrons oti C {< R) mostly 
ranges between 0.8 and 1.1. The distribution of secondary electrons, f es , is determined by the loss equation, which reads 

pesiP) = ~\fhs{p) b(p)] + Hp), (20) 

where p and b(p) have been previously defined and i(p) is the term describing the source of particles. For a steady state 
solution, if i(p) oc fr r and b(p) oc p* , then the particle distribution is 

/ es (p)cxp- r+1 - t . (21) 

We point out that most of the CR ions produced in our simulation have a spectrum characterized by a qi — 4.0 — 4.2 
(Miniati et al. 2001), which translates into a slope for the source function of secondary electrons close to the same value, 
i.e., r ~ 4.2. For IC losses, relevant for the HXR emitting electrons, t = 2; thus from eq. 13 

(r - 1 + t) - 3 

aic = ^ L1 ' t 22 ) 

in rough accord with the values in the top plots of Fig. 5. We notice, however, in those plots the presence of several points 
indicating values of «i C (< R) < 1. In addition, the number of such points becomes more numerous as we move toward 
the right panel, i.e., larger averaging volumes (Fig. 5). These values of the spectral index would imply the presence of 
secondary electrons with distributions, f es {p) oc p~ q "° ; flatter than q es = r — 1 + 1 = 4 — 1 + 2 = 5, i.e., than our simulated 
CR ions (with ~ 4) can generate. Such distributions are in fact produced by the re-acceleration that some secondary 
electrons undergo as they encounter a shock wave in the ICM. This event is more likely to happen at the outskirts of a 
cluster where shocks are more frequent, in accord with the trend of the top panels in Fig. 5. However, since we expect 
such re-accelerated secondary CRes to behave exactly as the primary electrons discussed above, the flattening of the 
spectral index should be only "visible" very close to the shocks and should disappear when the radiation from the large 
post-shock emission region is considered. Thus, due to the same coarse grid effects described above in the context of 
primary electrons, the spectral indexes reported in the top panels of Fig. 5 tend to be smaller than they should. Unlike 
the case of the primary electrons, however, correcting for this effect is highly impractical because it would require following 
(numerically) the re-accelerated particles as an independent population. In any case, we expect the values of a ic {< R) 
reported in the panel corresponding to an averaging volume of R = 0.3 Mpc to be only little affected by the re-accelerated 
populations of secondary CRes and, therefore, fairly representative of the simulation results. 

Finally we point out that for EUV emitting secondary electrons, with Lorentz factor 7 <~ 300, the loss mechanism are 
not only IC emission but also Coulomb collisions, implying a smaller value for both t and a. 

3.1.1. Hard X-ray Emission 

Having computed the fitting parameters summarized in Table 1 we can now estimate the HXR fluxes for Coma, Abell 
2199 and Abell 2256 clusters as predicted by our simulation model. For this exercise, an extrapolation has been applied 
using the parameters in Table 1, since the collapsed objects identified in simulation data have temperature < 3 keV while 
the observed clusters used for comparison have higher temperature. In order to define the scope of such extrapolations, 
first of all we notice that the the plots reported in Table 1 show quite a bit of (real) scatter, which we attribute to the 
different evolution of the collapsed objects. Thus our plots will only allow an estimate of average values and within a factor 
of the order of the scatter. As pointed out above, the scatter in the bottom-right panel of Fig. 3 is particularly large. With 
this clarification, in Table 2 we report the name of each observed cluster in the first column, and the measured red-shift, 
temperature, HXR flux and spectral index in columns 2-5 respectively (references in column 10). The expected average 
HXR flux and spectral index according to our secondary model are reported in columns 6 and 7 respectively. Finally, 
column 8 and 9 refer to the primary CRes model: the latter lists the spectral index associated with each population of 
CR particles and the former refers to the required R e / P value in order for the HXR flux from primary CRes to match the 
reported observations. 

The values reported in Table 2 indicate that the total number of secondary electrons is too small to account for the 
full measured HXR flux from any of Coma, Abell 2256 or Abell 2199. In fact, the calculated emission from secondary 
electrons in our simulation corresponds at most to a few percent of the reported measurements. Higher values of the 
total injection efficiency r/ in j (see §2.2) could allow a larger population of CR ions and, therefore, of secondary electrons. 
However, the resulting non-thermal pressure and 7-ray flux associated with the CR ions would then become inconsistent 
with the observational limits (Miniati et al. 2001; Blasi & Colafrancesco 1999; Blasi 1999). 

On the other hand, in order to account for reported detections in the HXR wave-band with the emission from the 
primary electrons in our simulation, would require R e / P > 0.2 for Coma and Abell 2256 and R e / P > 0.025 for Abell 2199. 
Since measurements for galactic CRs give R e / P ~ 0.01 — 0.05 (Mueller & Tang 1987), it would appear that only the case 
of Abell 2199 could be explained in terms of IC emission from primary electrons. However, as already pointed out, our 
plot in the bottom right panel of Fig. 3 exhibits a large scatter, of about one order of magnitude around the average 
values at each temperature, due to the different accretion/merger status of the simulated groups/cluster. Thus the HXR 
emission from Coma cluster and Abell 2256 is also consistent with the primary model, if a major accretion/merger process 
occurred in the recent history of these clusters, as indeed indicated by the observations (Briel et al. 1991; Burns et al. 
1994). 

As regarding the spatial distribution of the HXR radiation, it is unknown for the Coma cluster (Fusco-Femiano et al. 
1999), but there is some indication that it may be more widely distributed than the thermal emission for Abell 2199 
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(Kaastra et al. 1999), as it would result from the IC process (Sarazin & Lieu 1998). The emission from simulated 
secondary CRes is smoother and less extended than that from primary CRes, but the available information on the spatial 
distribution of detected HXR emission is not sufficient yet to discriminate observationally between the two alternative 
models. 

3.1.2. Extreme Ultra-violet Emission 

Regarding the EUV radiation, in Table 3 we compare the published EUV excess measurements in column 3 with our 
numerical results for: the EUV luminosity due to IC emission from secondary CRes in column 4; and the value of R e / P 
required to match the observations for primary CRes in column 5. 

The values computed in Table 3 for a cluster with the same temperature as Coma, indicate that both primary and, 
perhaps even, secondary CRes could produce through IC emission an EUV radiation flux comparable to the measurement 
reported by Bowyer et al. (1999). Our numerical model, however, falls short of the detected emission from Virgo cluster 
(Bcrghocfer et al. 2000; Lieu et al. 1996b) by about two orders of magnitude, although the presence of M81 may induce 
higher non-thermal particle populations in the ICM of Virgo and enhance the emission in this wave-band. For the 
remaining two clusters Abell 1795 and Abell 2199, our numerical model for both primary (assuming R e / p ~ 10~ 2 ) and 
secondary electrons, indicates a luminosity of order of a few 10 41 — 10 42 erg s _1 , well below the reported detections. The 
two order of magnitude difference is significant even considering the scatter in the reported plots and the limitations in our 
work. So for these two clusters we suspect that the reported EUV excesses are unlikely to be produced by IC scattering 
off cosmic microwave background photons of CRes shock-accelerated in the ICM. This conclusion is further supported 
by the following insight. The simulation results in Figs. 4 show a clear trend of the EUV luminosity increasing with 
the temperature of the collapsed object for both primary and secondary electrons. Such a trend is not present in the 
observational data. According to Table 3, the reported EUV luminosity for both Abell 1795 and Abell 2199 are larger by 
about two orders of magnitude than the value reported for Coma by Bowyer et al. (1999). This contradicts the theoretical 
trends shown in our plots and cannot be accounted for by the scatter among the simulated groups/clusters. We note 
that this peculiar feature of the observational results will affect any model in which the EUV radiation is produced by 
energetic electrons generated through general cluster evolution. In fact, whether accelerated by accretion shock waves, 
merger shock waves, intra-cluster turbulence or generated in hadronic interactions, the energy content of the relativistic 
electrons can always be expressed as a fraction of the total energy of the cluster (which scales with its temperature). 

The synthetic EUV images can also be compared to the thermal X-ray and synchrotron radio emission maps displayed 
on the same scale in Figs. 9 and 10. As we can see, the EUV emissivity from both primary and secondary electrons is more 
extended than the thermal X-rays, as observed in real clusters. However, the EUV radiation in these IC based models 
is also more extended than the radio emissivity; this contradicts observational results, as first pointed out by Bowyer & 
Berghoefcr (1998). Thus, the morphological properties of the emitting region pose another challenge to our attempt to 
attribute the EUV excess to CRes of cosmological origin. 

3.2. Synchrotron Radiation 

The synchrotron radiation emitted by CRes has been calculated using eq. 14. We remind the reader that we have 
assumed for this calculation a volume averaged magnetic field in a Coma-like cluster of order 3 /zG and an emitting 
volume of radius 1.3 /i _1 Mpc. Fig. 6 displays the simulated radio spectral power, P„, at 1.4 GHz (left) and 330 MHz 
(right) from secondary (top) and primary (bottom) electrons, respectively, as a function of the group/cluster temperature. 
For each plot, the best fit parameters for a power-law correlation (cf. §2.6) are reported in Table 1. Although the emitting 
electrons are approximately the same as those producing HXRs, the temperature dependence of P u is steeper now, 
reflecting the additional scaling of the magnetic field with the group/cluster temperature. The latter is close to T° 75 (see 
§3.2.1) in our simulation, and indeed the values of <j> have increased by approximately one unit with respect to the IC 
scattering emission. 

By analogy with the HXR emission, we have calculated the volume-averaged synchrotron spectral indices for radiation 
at 1.4 GHz and 31 MHz, which bracket the range of most commonly observed frequencies. These spectral indices are 
illustrated in Fig. 7 for secondary electrons and in Fig. 8 for the primary electrons. As for the HXR case, the spectral 
index associated to the emission from primary electrons represents the value averaged over the post-shock emitting volume. 
As discussed at length in §3.1, the spectral index of the CRes just behind the acceleration region should be smaller by one 
half with respect to what reported in those figures. First of all, those plots do not show any relevant difference between 
values of the spectral index at the frequencies of 1.4 GHz and 31 MHz, both in the case of primary or secondary electrons. 
This indicates that within the energy range that correspond to the above frequencies, there is no break due to radiative 
losses in the particle distribution. Similarly to and with analogous considerations as for the HXR IC emission, we find 
Oi S yn{< R) ^ 1-1 for the secondary electrons and 1 — 1.3 for the primary electrons in the large averaging volume which 
contains the acceleration sites. Again, this is because the particles responsible for the emission in the two cases (radio 
and HXR) represent approximately the same populations. 

3.2.1. Secondary Electrons and Radio Halos 

The secondary electrons from our numerical simulation can account for a number of features observed in radio halos 
including the morphology of the emitting region, the polarization properties and the temperature dependence of the 
synchrotron power at 1.4GHz (Liang et al. 2000). 
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In Fig. 9 we show synthetic maps of the surface brightness from synchrotron radiation from secondary electrons for 
the same object as in Fig. 1. In each of the three panels we present the radio emission at different wavelengths, namely 
1.4 GHz (top-left), 330 MHz (top-right), 31 MHz (bottom- left) . For reference, we also show the X-ray map produced 
by the thermal bremsstrahlung. The radio emission extends throughout the X-ray group/cluster, indeed resembling the 
spatial features of observed radio halos. By quantitative analysis of the emission map the lower frequency emission 
appears slightly more extended, again in accord with the observations. This trend is due to the decay of the magnetic 
field strength at large distances from the cluster centers and the fact that, according to cq. 15, the emissivity dims first 
at higher frequencies. A comparison with the images of X-ray from thermal bremsstrahlung reveals a higher level of 
substructure in the radio emission than in the thermal emission. Such sub-structure resembles to some extent that seen 
in the synthetic 7-ray images of Miniati et al. (2001) and is connected to the underlying shock distribution in the ICM. 
However, the synchrotron map will also depend on the structure of the magnetic field, which, in fully resolved simulations 
or in real clusters might extend to regions smaller than we have been able to simulate here. In any case, radio spatial 
features are also seen in real observations. There, the region responsible for the radio emission often appears irregular 
and, in this sense, shows some departure from the observed X-ray image (e.g., Abell 665, Giovannini & Feretti 2000; Abell 
520, Govoni ct al. 2001; and CI 0016+16, Giovannini & Feretti 2000). 

According to our calculations, a typical Coma-like cluster (T x — 8.3 keV) should exhibit a radio power at 1.4 GHz of 
order 1.5 x 10 24 W Hz -1 , close to the observations (see Table 4). We caution that this result depends on our choice of 
the magnetic field normalization of a few jiG for a Coma-like cluster. The significance of this finding is that the energy 
supplied by structure shocks would suffice to power radio halos if magnetic fields of a few fiG exist inside such clusters. 
In addition, our numerical findings indicate a scaling of Piaghz oc T 4 2 , in good agreement with the observational values 
(Liang et al. 2000). In fact, the synchrotron power in the secondary model should scale as 

r>l+a syn 

P syn cx iV es B 1+a °y" cx N cr n gas — ^ oc T 3 T 2 ^ oc T 3 + 2 ^ (23) 

B z /$ir + U cm b 

where a syn ~ 1 is the radio spectral index, N es cx N cr n gas /(B 2 /8ir + u cm b) is the number of secondary electrons, 
proportional to the product of the of cosmic ray ions and the gas density. The factor l/(B 2 /8n + u cm b) simply accounts 
in a steady state for the reduction of secondary electrons from radiative losses. Also, (\>b indicates parametrically the 
temperature dependence of the magnetic field strength, i.e., B cx T^ B . In §3.1.1 we already saw that N es cx N cr n gas cx T 3 . 
In addition, from our simulation we have found 4>b ~ 0.75. The amplification of the magnetic field is primarily due to 
shear flows developed during large scale structure formation and, to a lesser extent, to gas compression. The accuracy of 
the magnetic field strength dependence on the group/cluster temperature in the simulation is not clear, since the B — T x 
relation is observationally unknown, although now under investigation (Kronberg, priv. comm.). Clearly, though, the 
coarse grid effects limit the amount of shear and reduce the growth of the field, as it has already been pointed out in 
analogous studies (Rocttigcr ct al. 1999b; Dolag et al. 1999). 

With the above temperature scalings of N es and B, the result P\agh z oc T 4 2 is roughly retrieved. However, since 
otsyn ~ 1, when B > (87ru cm fc) 1 / 2 (l + z) 2 ~ 3.3y(iG, the synchrotron power becomes independent of the magnetic field 
strength and its dependence on the cluster temperature inevitably flattens. For the temperature scaling properties of the 
magnetic field in our simulation and the adopted normalization, such a turnover should occur at a temperature above 9.4 
keV (at z — 0). Given the poor statistics of real clusters with temperature beyond this value, conclusions in this respect 
are premature. 

We have also investigated the fractional polarization that characterizes the integrated synchrotron emission in the 
synthetic maps. From our investigation we find that the degree of polarization is very low throughout the simulated 
groups/clusters, typically of order of a few %. This result holds for the radiation at 1.4GHz and 330 MHz. In this 
calculation, no Faraday depolarization was included. Its inclusion would further reduce the degree of polarization found 
here (Dolag & Enfilin 2000), reinforcing our conclusions below. The low degree of polarization of synchrotron emission is 
one of the distinctive features of observed radio halos. The outcome in our model is probably due to the highly complex 
structure of the magnetic field inside the group/cluster. In fact, the radiation from different depths along the line of sight 
will imprint a polarization angle determined by the integrated structure of the field. Since the field is a varying function 
of space, the sum of all the contributions along one line of sight and orthogonal field projections of equal magnitude lead 
to a null net polarization, the emergent radiation from a disordered field will be mostly unpolarized. 

Another important quantity that is determined by the observations is the spectral index of the radiation. In the 
secondary CRe model this should be > 1.1. For the case of the Coma cluster, the best studied example, the integrated 
radio flux below 1.4 GHz is well fitted by a spectral index 1.16-1.3 (Deiss et al. 1997). This is marginally consistent with 
our results and, in fact, our current model underestimates the radio emission of the Coma cluster at 330 MHz with respect 
to the data of Venturi et al. (1990). Additionally, the observed spectral index shows a steepening with increasing radius. 
No conclusion can be drawn from our plots in this respect. In fact, although Fig. 7 shows that the emission averaged 
spectral index tends to decrease slightly when larger volumes are considered, this is a spurious effect as was discussed in 
§3.1. In general the observational values of the spectral index may be slightly higher than what we predict (Hanish 1982). 
In the context of secondary electrons from cosmic shocks this would correspond to slightly weaker shocks in the ICM than 
those typically found within group/cluster ICMs in our simulation. In fact, there is a possible tendency for the Mach 
numbers we compute to be slightly overestimated, especially in regions close to the cluster core where with moderate 
spatial resolution the shock compression can be "confused" by the gravitational compression (see Miniati et al. 2000, 
for details). Although a precise estimate should be carried out through higher resolution simulations, a simple estimate 
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indicates that 

^-^-"ir^^Bf^ < 24 > 

where "f gas = 5/3 is the adiabatic index of the gas. Therefore, an error AM ~ 1 in the Mach number estimate translates 
into an error Aa ~ —0.2 for the spectral index, which would allow for spectra steeper than reported in our plots. 

3.2.2. Primary Electrons 

In Fig. 10 we present synthetic maps of the surface brightness due to synchrotron emission from primary electrons at 
1.4 GHz (top-left), 330 MHz (top-right) and 31 MHz (bottom-left) for the same object of Figs. 9. On the bottom-right 
panel, we also present on the same scale a two-dimensional cut of the contour levels of shock compression in a plane 
through the center of the collapsed object and perpendicular to the line of sight of the other three panels. In contrast to 
the emission produced by the secondary electrons, the primary electron emission maps are strikingly inhomogeneous, with 
the primary electrons (and emission) strongly concentrated around shocks. Fig. 10 resembles the HXR images presented 
in Fig. 1 (right panel), which also refers to the same object. 

The location of these synchrotron emission regions at the immediate outskirts of the group/cluster together with their 
irregular and extended morphology, may suggest a connection with the observed radio relics. The suggestion is reinforced 
by the finding that the fractional polarization through the extended regions is quite high, of order 30-50 %, comparable 
to the values suggested by the observations. Unlike the the secondary electrons, the primary electrons produce emission 
primarily in the immediate post-shock regions. There the magnetic field is highly aligned due to the enhancement of 
the perpendicular component by the shock compression. Thus the radiation suffers only minor structural depolarization 
along the line of sight and emerges highly polarized. We do point out that our results might overestimate the polarization 
expected in real clusters because finite resolution limits the possible complexity of the magnetic field structure. 

In Table 5 we present a comparison between the observed properties of radio relics for a few sources and the averaged 
values predicted in our numerical work. To define the scope of this comparison, we point out again that the predictions 
of our model, illustrated in the bottom panels of Fig. 6, indicate a range of values which span a factor of a few about 
some average. With this in mind, we notice that in all of the presented cases the mean predicted total power at 1.4 GHz 
is consistent with the observed values for a conservative estimate of R e / P of 10 -2 — 10~ 3 . This result signifies that, again, 
the power supplied at shocks by accretion flows is sufficient to generate the emission observed in extended radio relics, for 
conservative values of shock acceleration injection and efficiency. 

We have also calculated the volume-averaged spectral index (cf. §3.1) for each pixel in the images of Figs. 10. We find 
that in the brightest regions a syn ~ 1 — 1.3 at both 1.4 GHz and 31 MHz. Where the emissivity is low, the spectral index 
steepens to values a syn « 1.4 — 1.5 (at 31 GHz). These values are typical for all the groups/clusters analyzed with the 
simulation data (Miniati 2000) and are indeed consistent with the observed values for relics, a few of which are reported in 
table 5. In addition, the above volume-averaged values of a syn corresponds to a spectral index of order ~ 0.5 — 0.8 for the 
bright emission regions associated with the freshly accelerated populations just behind the shock (cf. 3.1). Observation 
of the relic source 2006-56 is consistent with the above depiction as the spectral index increases from a syn si 0.5 from 
the edge of the source to about a syn ~ 1.1 toward its center (Rottgering ct al. 1997). This trend may indicate that the 
emitting particles in that source are indeed accelerated by shocks, as was pointed out by Enfilin et al. (1998). 



4. DISCUSSION AND CONCLUSIONS 

In this paper we have explored the acceleration, losses and spatial transport of cosmic ray electrons for the first time 
in a fully cosmological simulation of large scale structure formation. In our simulation, CRes were either directly injected 
at shocks (primary) or generated in hadronic collisions of CR ions and the nuclei of the ICM (secondary) . In general the 
populations of CRs thus generated are described by power laws, f(p) cx p~ q , with index q depending on the strength of 
the shock. Our results indicate that the total number of relativistic electrons present in each group/cluster scales with 
the temperature of the ICM in the sense that more numerous populations are present in hotter groups / clusters. 

Non-thermal emissions due to relativistic electrons have been computed at different wavelengths in order to address 
the issues related to the recent EUV/HXR observations of clusters of galaxies. As for the excess of emission in the 
EUV wave-band, we find that the electrons generated in our simulation arc probably insufficient to produce the claimed 
detections (Lieu et al. 1996b, a; Mittaz et al. 1998; Berghoefer et al. 2000), except for the case of Coma cluster reported 
by Bowyer et al. (1999). Furthermore, the EUV excess reported so far has shown no trend with cluster temperature, 
contrary to the results of our simulation. Absence of such trends would suggest that processes responsible for generation 
of non-thermal electrons do not depend on the energetics of groups/clusters. 

The contribution to the HXR emissivity from the secondary electrons amounts to only a few percent of the currently 
reported values (namely, for the Coma cluster and Abell 2556 in the bandwidth 20-80 keV and Abell 2199 in the bandwidth 
0.13-100 keV). Increasing the total number of secondary electrons enough to account for all of the HXR excess emission 
would imply a corresponding increase of parent CR ions and a consequent 7-ray flux much above the upper limits 
established by EGRET (Blasi & Colafrancesco 1999; Blasi 1999; Miniati et al. 2001). 

The average HXR emission from primary electron extrapolated from our plot in the bottom right panel of Fig. 3 
indicates a contribution from this component of about 10% or so of the reported detections in this wave-band. This plot, 
however, is characterized by a large spread of about an order of magnitude around the average, a reflection of the different 
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level of shock activity in the ICM of individual groups/clusters. Thus, the reported HXR emission from Coma cluster 
and Abell 2556 could be accounted for by IC emission of primary electrons, given that these clusters are known to be 
respectively in a post-merger phase (Burns et al. 1994) or in the process of merging (Bricl ct al. 1991). 

From the analysis of the synchrotron emission, we find that our secondary CRes could explain several features of radio 
halos; namely the total power at 1.4 GHz, the fractional polarization, and the morphology and, marginally, even the 
spectral index. The radio power, Piaghz, as reported by Liang et al. (2000), for clusters at different temperatures is 
produced by the secondary electron model, provided the magnetic field varies with the group/cluster temperature as 
steeply as found in our simulation. Consistency here is achieved by assuming a normalization for the magnetic field that 
gives the volume- averaged value of 3 /xG for a Coma- like cluster, as inferred from Faraday rotation measure estimates 
(Kim et al. 1989, 1990). Numerous properties of radio halos are reproduced quite naturally according to our simulation, 
although the spectral index reported by observations is sometimes steeper than found here (see also Chi et al. 1996). 
Future higher resolution simulations should allow us to determine whether or not this is a numerical artifact, as suggested 
in §3.2.1. A similar result about the secondary model for radio halos was recently obtained by Dolag & Enfilin (2000), 
although in their case the CR population was not simulated but assumed from the beginning. In our simulation the 
spectra of accelerated particles are determined consistently with the cosmic shocks responsible directly or indirectly for 
their acceleration. Thus, particle distribution and spectra are not free parameters of the computation. 

Finally, we have shown that the primary electrons directly accelerated at cosmic shocks generate radio maps that 
are similar in radiated power, morphology and polarization fraction to observed radio relics. The accelerating shocks 
have been shown to be strong enough to produce quite flat populations of CRes in the region immediately adjacent to 
the acceleration sites. Our limited numerical resolution prevents us from reproducing correctly the spatial distribution 
of electron populations of different "ages" propagating away from the shock. However, the computed flat spectra of the 
shock-accelerated electron distributions indicate that the volume-averaged spectral index of the whole post-shock emission 
should be in the observed range, a = 1 — 1.5. 

Thus the energetics of accretion/merger shocks are in principle sufficient to power both the radio-relic emission and the 
reported detections of HXR excess. However, in order for the same electron population to account for the emission in 
both wave-bands for either Coma cluster or Abell 2256 (see table 2 and 5), it is required a volume-averaged magnetic field 
strength about an order of magnitude smaller than what assumed here. It is conceivable that the magnetic field strength 
drops by this amount as we move from the cluster center to its outskirts where the accretion shocks powering both the 
radio-relics and the HXR emission, would be located according to our simulations. And in fact this is supported by the 
magnetic field maps of our simulations (Miniati 2000). However, due to the limited numerical resolution such properties 
of the magnetic field cannot be considered as conclusive. 

It is important to mention the existence of alternative models for the existence of radio relics. In particular, it has been 
pointed out that radio relics could be sites where fossil radio plasma, expelled by radio galaxies at some point during 
cosmic history, is currently being revived through the encounter with a shock wave by means of adiabatic compression 
(Enfilin 1999; Enfilin & Gopal-Krisna 2001). The idea is compelling particularly because the ejection of relativistic plasma 
by radio galaxies in the intergalactic medium is an observational fact. It is shown that radio plasma up to 2 Gyr of age can 
be revived by strong shocks (Enfilin & Gopal-Krisna 2001). The statistics, i.e., the probability of actually seeing revived 
radio plasma at a given epoch, is not an easy task because it depends on the luminosity function of radio galaxies as well 
as the subsequent evolution of the radio plasma in cosmic environment. However, future low frequency observation should 
improve substantially our knowledge of the distribution of old radio plasma in cosmic environment (Enfilin & Gopal-Krisna 
2001). Nevertheless, there exist some important differences between radio relics produced by primary electrons accelerated 
at shocks and those resulting from revived radio plasma. In fact, according to recent magetohydrodynamical simulations 
of the interaction of radio plasma with a shock wave, a cloud of radio plasma becomes unstable upon being shocked, and 
develops into a complex filamentary structure (Enfilin & Bruggen 2001). It is not clear that such morphological properties 
belong to giant radio relics observed at the outskirts of cluster of galaxies. However, recent radio observations have shown 
that a sample of radio relics with sizes up to a few hundred kpc observed at very high resolution in a number of galaxy 
clusters possess just these morphological features (Slee et al. 2001). It is possible then that both mechanisms, namely 
shock acceleration and compression of old radio plasma, are at work in the ICM producing a large variety radio sources 
with distinct properties, as outlined above. This possibility, just speculated here, should be investigated in the future. 

From a broader perspective, it is very important to understand the actual origins of radio halos and relics, so one can 
assess the level of non-thermal activity in the ICM. If, as provided in the present model, radio halos are produced by 
secondary electrons, they reveal the existence of a related population of CR ions that have developed during the formation 
of the large scale structure (Volk et al. 1996; Berezinsky et al. 1997; Miniati 2000; Miniati et al. 2001). Such populations 
may contain a significant fraction of the total energy density of the ICM and could be, therefore, dynamically important 
(Miniati 2000; Miniati et al. 2001), with many cosmological implications. 
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Table 1 

Temperature scaling relations 





secondary e 




primary a e 




Energy Band 


K 




K 






(erg s" 1 cm -2 ) 




(erg s _1 cm -2 ) 




F ic HXR: 20-80 keV 


2.8 x 10" 13 


2.9 


8.0 x 10~ n 


1.5 


F ic HXR: 0.13-100 keV 


7.0 x 10~ 12 


2.9 


3.6 x 10- 9 


1.9 




(erg s _i ) 




(erg s _i ) 




L ic EUV: 65-248 eV 


2.2 x 10 41 


2.9 


3.8 x 10 43 


1.9 




(W Hz" 1 ) 




(W Hz^ 1 ) 




f 1.4GHz 


6.3 x 10 23 


4.2 


7.3 x 10 25 


2.6 


-P330MHz 


2.4 x 10 24 


4.1 


2.1 x 10 26 


2.6 


-P74MHZ 


1.2 x 10 25 


4.2 


7.4 x 10 26 


2.7 


-P3IMHZ 


2.9 x 10 25 


4.2 


1.6 x TO 27 


2.8 



a The values reported in the table for primary electrons assume R e / P = 1. 



Table 2 

Flux of HXR radiation: measurements vs predictions. 











predictions 








measurements 




secondaries primaries' 1 




Cluster 


z 


T x h F H xr c 


o f 


Fhxr c a lc f i? e/p g a ic h 


reference (s) 


Coma d 


0.0232 


8.3 2.2 


0.7-3.6 


5.2 x 10~ 2 -1.05 0.2 1-1.3 


Fusco-Fcmiano et al. (1999) 


Abell 2256 d 


0.0581 


6.95 1.2 


0.3-1.7 


4.8 x 10~ 3 -1.05 0.8 1-1.3 


Fusco-Femiano et al. (2000) 


Abell 2199 c 


0.0305 


4.6 3.2 


1.5-2.0 


0.13 -1.05 0.03 1-1.3 


Kaastra et al. (1999) 



"We report here the required value for R e / P in order to explain the observed emissivity with the primary electrons model. 

b Temperatures are measured in keV. 

c Fluxes are measured in units 10 — 11 erg s — 1 cm" 2 . 

d For Coma cluster and Abell 2256 all fluxes correspond to the bandwidth 20-80 keV. 
e For Abell 2199 cluster all fluxes correspond to the bandwidth 0.13-100 keV. 
■^The spectral index cti c is defined in §13. 

9 We report here the required value for R e / p in order to explain the observed emissivity with the primary electrons model. 

h The reported values of the spectral index cti c refer to averages over the post-shock emitting volume (cf. §3.1). For freshly accelerated particles 
just behind the shock they should be in the range 0.5 — 0.8. 
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Table 3 

Emission of EUV radiation in the EUV band 65 - 248 eV: measurements vs predictions. 

Predictions' 1 



Cluster 


T x 


measurements 


secondaries 




reference(s) 




keV 


(erg s" 1 ) 


(erg s" 1 ) 




Coma (A 1656) 


8.3 


1.5 x 10 42 


4.0 x 10 41 


6 x 10" 3 


Bowyer et al. (1999) 




8.3 


5.0 x 10 43 


4.0 x 10 41 


0.2 


Lieu et al. (1996b) 


Virgo 


1.8 


5.2 x 10 42 


4.8 x 10 39 


0.4 


Berghoefer et al. (2000) 




1.8 


9.0 x 10 42 


4.8 x 10 39 


0.7 


Lieu et al. (1996b) 


Abell 1795 


5.8 


n. d. 


1.4 x 10 41 




Bowyer et al. (1999) 




5.8 


2.0 x 10 45 


1.4 x 10 41 
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Mittaz et al. (1998) 


Abcll 2199 


4.5 


n. d. 


7.0 x 10 40 




Bowyer et al. (1999) 




4.5 


5.0 x 10 43 


7.0 x 10 40 


0.7 


Sarazin & Kempner (2000) 



a We report here the required value for R e / p in order to explain the observed emissivity with the 
primary electrons model. 



Table 4 

Radio power at 1.4 GHz: measurements vs predictions of secondary electrons model. 



Quantity 


measurements 


model 


reference (s) 


Pi.4GH,(8.3keV) a 


1.5 


1.5 


Deiss et al. (1997) 


4> 


4-5.2 


- 4.2 


Liang et al. (2000) 




0.5-< 1.6 


>1-1 


Hanish (1982); Deiss et al. (1997) 


Pol. 


< 10 % 


a few % 


Feretti (1999) 



a The radio power is measured in units 10 24 W Hz 1 . Here we compare the normalizations of the P\AGHz vs T x relation by taking the measured 
and predicted values of PiAGHz f° r the Coma cluster. 
'The range includes values for different clusters as well as spatial variations inside the same cluster. 



Table 5 

Properties of Radio Relics: measurements vs predictions of primary electrons model. 











measurement 


modcl b 




Relic 


Cluster 


freq. 


Flux a 


Pol. 




a 


p c 


ref.(s) 


1253-275 


Coma 


1.4 GHz 


0.16 


25-30 % 


1. 


1 ± 0.2 d 


1 x 10~ 3 


Venturi et al. (1990); Giovannini et al. (1991) 






330 MHz 


1.4 


n. a. 


1. 


1±0.2 


2 x 10~ 3 


Giovannini et al. (1991) 


2006-56 


A 3667 


1.4 GHz 


0.7 


n. a. 


0. 


5 - 1.2 


4 x 10~ 2 


Goss et al. (1982); Rottgering et al. (1997) 


Comp. G 


A 2256 


1.4 GHz 


0.15 


20 % 


0. 


8 ±0.1 


6 x 10" 3 


Rottgcring et al. (1994) 


Comp. H 


A 2256 


1.4 GHz 


0.1 


20 % 


0. 


8 ±0.1 


5 x 10" 3 


Rottgering et al. (1994) 


1140+203 


A 1367 


1.4 GHz 


0.18 


n. a. 


1. 


9 ±0.5 


5 x 10~ 3 


Gavazzi & Trinchieri (1983) 


The radio flux 


is measured 


in Jy = 10 - ' 


!3 erg cm 


-2 











''The spectral index in the model is always in the range 0.5-0.8 for the emitting region immediately behind the shock and 1-1.3 when the 
volume-average value is considered (see §3.1 for further details). Likewise, the prediction for the fractional polarization is always 30-50 %, 
although geometrical effects may reduce the observed value with respect to it (Enfilin et al. 1998). 

c Here we report the required value for -R e / p in order to explain the observed emissivity with the primary electrons model. 

d Andcrnach ct al. (1984) find 1.18. 
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Fig. 1. — Emission map of HXR in units 2 X 10 4 erg cm 2 s 1 per pixel, from secondary (left) and primary (right) electrons at 50 keV 
for one of the collapsed objects generated in the simulation. The linear size of the panels correspond to 5.5 h -1 Mpc. 
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Fig. 2. — Emission map in units 4 X 10—7 erg cm 2 s 1 per pixel, of EUV radiation for the same object as in Fig. 1. Left and right panels 
correspond to secondary and primary electrons, respectively. The linear size of the panels correspond to 5.5 h -1 Mpc. 
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Fig. 3. — Group/cluster HXR flux of IC emission in 0.13-100 keV (left) and in 20-80 keV (right) due to secondary (top) and primary 
electrons (bottom), respectively, as a function of the group/cluster core temperature. The emissivity is integrated over a spherical volume 
within a radius of 1.3 /i _1 Mpc. Note that the physical values of the luminosities for primary electrons should be scaled with the electron to 
proton injection ratio, R e / p . 
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Fig. 4. — Group/cluster EUV luminosity from IC emission due to secondary electrons (left) and primary electrons (right) as a function 
of the group/cluster core temperature. The emissivity is integrated over a spherical volume within a radius of 1.3 /i _1 Mpc. Note that the 
physical values of the luminosities for primary electrons should be scaled with the electron to proton injection ratio, R e / p . 
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Fig. 5. — Top: Volume-averaged spectral index for the HXR radiation at 50 keV from secondary electrons as a function of the group/cluster 
core temperature for four different averaging volume radii, i.e., from left to right R = 0.3,0.5,0.9, 1.3 7i -1 Mpc. Bottom: same as top but now 
for primary electrons. 
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Fig. 6. — Group/cluster synchrotron power at 1.4 GHz (top-left) and at 330MHz (top-right) from secondary electrons and at 1.4 GHz 
(bottom-left) and at 330MHz (bottom-right) from primary electrons as a function of the core temperature. The emissivity is integrated over 
a spherical volume within a radius of 1.3 7i -1 Mpc. Note that the physical values of the luminosities for primary electrons should be scaled 
with the electron to proton injection ratio, R e /p. 
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Fig. 7. — Volumc-avcragcd spectral index for the synchrotron radiation from secondary electrons at 1.4 GHz (top) and 31 MHz (bottom) 
as a function of the group/cluster core temperature. The four panels in each row correspond to the following four different averaging volume 
radii: from left to right R = 0.3, 0.5, 0.9, 1.3 /i _1 Mpc. 
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Fig. 8. — As in Fig. 7 but now for primary electrons. 
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Fig. 9. — Radio map of surface brightness in units of Jansky per pixel (1 Jy = 10 -23 erg cm -2 ) for the same object as in Fig. 1, at 1.4 
GHz (top-left), 330 MHz (top-right) and 31 MHz (bottom- left). The radiation is synchrotron emission generated by secondary electrons in 
the group/cluster magnetic field. For comparison we also show the thermal X-ray map (bottom-right) in units erg cm~ 2 s _1 . The linear size 
of the panels correspond to 5.5 /i -1 Mpc. 
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Fig. 10. — Same as for Fig. 9 but for primary electrons in the top-left, top-right and bottom left panels. A two-dimensional cut of 
the contour levels of shock compression in a plane through the center of the collapsed object is also shown on the bottom-right panel for 
comparison. 
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